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ABSTRACT 



We present the results of searches for high-energy muon neutrinos from 41 gamma- 
ray bursts (GRBs) in the northern sky with the IceCube detector in its 22-string con- 
figuration active in 2007/2008. The searches cover both the prompt and a possible 
precursor emission as well as a model-independent, wide time window of — lh to +3h 
around each GRB. In contrast to previous searches with a large GRB population, we do 
not utilize a standard Waxman-B ahead GRB flux for the prompt emission but calcu- 
late individual neutrino spectra for all 41 GRBs from the burst parameters measured by 
satellites. For all three time windows the best estimate for the number of signal events 
is zero. Therefore, we place 90% CL upper limits on the fluence from the prompt phase 
of 3.7 x 10~ 3 erg cm -2 (72TeV - 6.5 PeV) and on the fluence from the precursor phase 
of 2.3 x 10~ 3 ergcm~ 2 (2.2 TeV - 55TeV), where the quoted energy ranges contain 90% 
of the expected signal events in the detector. The 90% CL upper limit for the wide time 
window is 2.7 x 10 _3 ergcm~ 2 (3 TeV - 2.8 PeV) assuming an E~~ 2 flux. 

Subject headings: gamma-ray bursts: general - methods: data analysis - neutrinos - 
telescopes 



1. Introduction 

Gamma-ray bursts (GRBs) are among the most violent events in the universe and among 
the few plausible candidates for sources of the ultra-high energy cosmic rays. So-called long- 
duration GRBs (> 2 s) are thought to originate from the collapse of a massive star into a black hole 
(Woosley 1993), whereas short-duration GRBs (< 2 s) are believed to be the result of the merger 
of two compact objects (e.g., neutron stars) into a black hole (Eichler et al. 1989). Though quite 
different in nature both scenarios are consistent with the currently leading model for GRBs, the 
fireball model (Meszaros & Rees 1993), with the energy source (central engine) being the rapid 
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accretion of a large mass onto the newly formed black hole. In this model, a highly relativistic 
outflow (fireball) dissipates its energy via synchrotron or inverse Compton radiation of electrons 
accelerated in internal shock fronts (Narayan et al. 1992; Rees & Meszaros 1994; Sari & Piran 
1997). This radiation in the keV-MeV range is observed as the 7-ray signal. In case of long GRBs 
the energy in gamma rays is typically of (9(10 51 -10 54 erg x 0,/Att) where Q is the opening angle for 
the 7-ray emission. Short GRBs are observed to release about a factor 100 less energy. 

In addition to electrons, protons are thought to be accelerated via the Fermi mechanism, 
resulting in an E~ 2 power law spectrum with energies up to 10 20 eV (Waxman 1995; Vietri 1995). 
The normalization of the proton spectrum is usually given in relation to the energy in electrons. 
The latter is linked to the energy in 7-ray photons through the synchrotron and inverse Compton 
energy-loss mechanisms. Protons of O(10 15 eV) interact with the keV-MeV photons forming a A + 
resonance which decays into pions (Waxman & Bahcall 1997). In the decay of the charged pions, 
neutrinos of energy 0(1O 14 eV) are produced with the approximate ratios (v e :v^:v T ) = (1:2:0) 1 , 
changing to about (1:1:1) at Earth due to oscillations (Learned k, Pakvasa 1995; Athar et al. 2006). 
First calculations of this prompt neutrino flux (Waxman &: Bahcall 1997; Alvarez-Muniz &i Halzen 
1999) used average GRB parameters and the GRB rate measured by BATSE to determine an 
all-sky neutrino flux from the GRB population. The AMANDA-II neutrino telescope (Achterberg 
et al. 2007, 2008) performed searches for this so-called Waxman-Bahcall GRB flux or similar GRB 
fluxes (Achterberg et al. 2007, 2008) with negative results. 

In a similar way, so-called precursor neutrinos can be generated when the expanding fireball is 
still inside the progenitor star (Razzaque et al. 2003). In this case, the accelerated protons interact 
with matter of the progenitor star or synchrotron photons. However, due to the large optical depth 
the synchrotron photons cannot escape the fireball and, hence, no 7-ray signal is observed. The 
time delay between the start of this neutrino emission and the prompt 7-ray signal is expected to 
be about 100 s. 

Observations of the early and late afterglow phases reveal that a large fraction of GRBs show 
X-ray flares superposed on the decaying light curve. Sometimes these flares are interpreted as a 
restart of the central engine that already generated the prompt emission (Burrows et al. 2007). If 
this is true, neutrino production with a similar spectrum as the prompt emission can be expected 
in the afterglow phase up to 10 4 s after the 7-ray signal (Murase & Nagataki 2006). Furthermore, 
production of neutrinos with energies around 10 18 GeV is expected when the shock fronts collide 
with the interstellar medium or the progenitor wind (Waxman & Bahcall 2000). 

In our analysis we search for muon neutrinos from GRBs recorded by satellites between 2007 
June 1 and 2008 April 4 in all three phases. For the prompt phase we utilize both an unbinned 
likelihood and a binned method. We find that the unbinned likelihood method has a significantly 
better discovery potential and is therefore used to obtain the limits presented in this paper. For 



1 Here and throughout the rest of the paper v denotes both neutrinos and antineutrinos. 



-6- 



searches in other emission phases we perform only an unbinned search. The paper is structured 
as follows: In Section 2 we define the neutrino spectra used for the different phases, followed by a 
description of the IceCube detector in Section 3. Afterwards, in Sections 4 and 5 the data sets and 
simulations are discussed, respectively In Section 6 the unbinned likelihood and binned methods 
are described and their performance is compared. Section 7 then presents the results followed by 
a discussion of systematic uncertainties. Finally, Section 8 sets the results into context with other 
observations. 

2. GRB neutrino-spectra and time windows 

The searches in this paper rely on the directional, temporal and spectral information obtained 
from satellite-based 7-ray observations which are distributed via the Gamma-ray burst Coordinate 
Network (GCN, NASA (2009)). Primarily, this information comes from Swift (Burrows et al. 2005) 
(also X-ray and UV observations), but also from Konus-Wind (NASA 1994), Super AGILE (Tavani 
et al. 2008), Integral (Mereghetti 2004), and other satellites of the Third Interplanetary Network 
(IPN 1990). In our analyses we only consider bursts which occurred at a declination above —5°. 
The southern sky is dominated by downgoing muons created by cosmic ray interactions with the 
atmosphere. By restricting our searches to the northern sky, the background from downward going 
muons is drastically reduced. The resolution of the GRB position from the satellites is better than 
0.1°, well below the resolution of the IceCube detector. It is therefore neglected in these analyses. 

2.1. Prompt emission 

We calculate the expected prompt neutrino spectrum in the internal shock scenario of the 
fireball model following Guetta et al. (2004) which is based on Waxman k, Bahcall (1997). For 
reference we list all formulae used in our calculations in Appendix A together with a definition of 
the various parameters. Prompt neutrino emission from GRBs is the result of meson production 
in collisions of accelerated protons and the observed 7-rays in the keV-MeV range. It is therefore 
expected to occur during the same time frame as the 7-ray emission and to track the photon 
energy spectrum. This is reflected in the similar functional form of F 7 (Equation (Al)) and F v 
(Equation (A3)). Due to the A + resonance condition, the neutrino energy is predicted to be 
inversely proportional to the photon energy, which is illustrated in the definition of a v and /3„ in 
Equation (A6). The further break in the neutrino spectrum above an energy €2 is due to synchrotron 
cooling of high energy pions and muons before producing neutrinos. The expected energy fluence in 
neutrinos is directly proportional to the measured energy fluence in photons (Equation (A8)). Here, 
some of the measured photon indices lead to diverging 7-ray spectral integrals if integrated from 
zero to infinity in energy. As the photon spectrum will not follow a broken power law spectrum to 
arbitrarily high or low energies, we limit the integration range for all GRBs from 1 keV to 10 MeV 
for which broken power law spectra have been observed by 7-ray satellites. 



Table 1. Burst Parameters, 7-Ray and Neutrino Spectra of All 41 GRBs (for Definitions of Parameters See Appendix A) 
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Note. — Columns: Tq - trigger time of satellite [UT], RA - right ascention of GRB [°], Dec - declination of GRB [°], T\ — Tq - start of prompt window [s], T2 — Tq - end 
of prompt window [s] , / T [MeV -1 cm -2 ], e-y [MeV] , f u [GcV - 1 cm -2 ], e u \ [PcV] , 6^2 [PcV] . The parameters f~ f and f v arc the fluxes at e~ ( and e v ^\ of the gamma-ray and 
neutrino spectrum, respectively (sec also Appendix A). 



* Parameter has not been measured. Instead, an average value is used (see Table 2). 
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Table 2. Average Values of GRB Parameters Taken from Becker (2008) 



Parameter 


Average Value 


A 


1.3MeV _1 cm- 2 a 


z 


2 




0.2 MeV 




1 


/3 7 


a 7 + 1 


t iso* 


10 51 erg cm -2 


1 jet 


300 


t * 

^var 


0.01s 


f * 


0.1 




0.1 


fe* 


0.1 



Note. - For parameter definitions see Ap- 
pendix. 

Corresponds to T 1 = 10~ 5 erg cm -2 between 
lOkeV and 10 MeV. 

*Not measured for any GRB. 



10 3 1 4 1 5 1 $ 10 7 1 8 1 ! 
Ev [GeV] 



Fig. 1. — Calculated neutrino spectra for all 41 GRBs (thin solid lines) compared to the standard 
Waxman-Bahcall spectrum for a single burst (thick dotted line). Also shown are the sum of all 41 
individual spectra (thick solid line) and the sum of 41 Waxman-Bahcall-like spectra (thick dashed 
line) . 

The time window of the prompt emission is determined from the information published in 
the GCN circulars and reports and checked against the measured 7-ray emission curves (available 
from NASA (2009)). In case of early or late emission outside the window specified by the satellite 
experiments the window is extended accordingly. The exact window definitions for the bursts are 
listed in Table 1. In previous publications on stacked searches for neutrinos from GRBs (Achterberg 
et al. 2007, 2008) it was assumed that the sum of all GRB spectra follows a Waxman-Bahcall GRB 
spectrum (Waxman & Bahcall 1997). However, the burst parameters can vary significantly from 
burst to burst (Guetta et al. 2004; Becker et al. 2006) and the GRB population used here (mostly 
observed by the Swift satellite) is different from the BATSE population on which the calculations 
from Waxman and Bahcall were based. Therefore, we use the measured parameters to calculate the 
neutrino spectra for each GRB individually as has already been done in the case of single, bright 
bursts (Stamatikos et al. 2005; Abbasi et al. 2009b). The parameters (see appendix A) for each 
GRB are listed in Table 1. In case a parameter has not been measured for a particular burst the 
average value listed in Table 2 is used. The table also contains assumed values of parameters that 
have not been measured for any of the GRBs. The resulting neutrino spectra for all 41 GRBs are 
shown in Figure 1 together with the standard Waxman-Bahcall spectrum 2 . Clear differences in the 



2 The fluence of a standard Waxman-Bahcall burst is calculated from the flux quoted in Waxman (2003) (Equation 
(17), given as a "diffuse" all-sky flux in GeV -1 sr _1 s _1 ) by multiplying it with 47rsr x f yr and dividing it by the 
assumed number of bursts per year (667). The resulting fluence is divided by two to account for neutrino oscillations 
(full mixing assumed). 
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shapes of the two summed spectra are observed together with an O(10) times lower overall fluence 
for the individual spectra. The latter is caused by the much higher sensitivity of Swift compared 
to BATSE. The observed differences stress the importance of using individual fluences in analyses. 

2.2. Precursor emission 

In the case of precursor emission, the neutrino-producing interactions occur while the fireball 
is still opaque to electromagnetic emission. Therefore, no analogous photons are observed and 
modeling of the emission on a per burst basis is not possible. We therefore use the fluence derived 
by Razzaque et al. (2003) for H-progenitor stars and assume that all 41 GRBs have such a precursor 
phase with the same fluence. The spectrum is significantly softer than that of the prompt emission. 
Below lOTeV, it follows an E~ 2 power law spectrum and has a sharp drop around 60TeV (see 
Figure 8). Neutrinos below 60TeV are mainly produced in interactions of accelerated protons with 
cold stellar protons whereas those above 60TeV originate from proton interactions with photons 
in the jet. The time window is taken as the 100 s immediately preceding the prompt time window. 
The window is chosen to be large enough to encompass the predicted emission and a potential delay 
between the phases. 

2.3. Wide window emission 

While specific predictions have been made for neutrino emission both before, during, and 
after the observed 7-ray emission of GRBs, there are many unknown quantities that factor into 
the calculation of fluence. It is therefore important to search for generic emission of high energy 
neutrinos in a reasonable time window surrounding the observed bursts. For the first time, we 
perform such a search in a wide time window (— lh to +3h) around each burst. The size of the 
window is motivated by possible precursor and afterglow emission, and limited by the requirement 
to keep backgrounds low. Rather than attempt to model the emission, we assume a generic E~ 2 
energy spectrum. Such a spectrum is in agreement with the assumed parent cosmic-ray spectrum 
and is distinguishable from the atmospheric neutrino background (see Section 5). 

3. Detector and data acquisition 

IceCube (Achterberg et al. 2006), the successor of the AMANDA experiment and the first next- 
generation neutrino telescope, is currently being installed in the deep ice at the geographic South 
Pole. Its final configuration will instrument a volume of about 1 km 3 of clear ice in depths between 
1450 m and 2450 m. Neutrinos are reconstructed by detecting the Cherenkov light from charged 
secondary particles, which are produced in interactions of the neutrinos with the nuclei in the ice 
or the bedrock below. The optical sensors, known as Digital Optical Modules (DOMs), consist 
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of a 25 cm Hamamatsu photomultiplier tube (PMT) housed in a pressure-resistant glass sphere 
and associated electronics (Abbasi et al. 2009a). They are mounted on vertical strings where each 
string carries 60 DOMs. The final detector will contain 86 such strings spaced horizontally at 
approximately 125 m intervals 3 . Physics data taking with IceCube started in 2006 with 9 strings 
installed. The completion of the detector construction is planned for the year 2011. The analyses 
described here use data taken with the 22-string configuration of the detector, which operated 
between 2007 May 31 and 2008 April 5. 

The data acquisition (DAQ) system of IceCube (Abbasi et al. 2009a) is based on local coinci- 
dences of photon signals (hits, threshold 0.25 photo-electrons) in neighboring or next-to-neighboring 
DOMs on a string within 1 fis. All data from DOMs belonging to a local coincidence are read out 
and the digitized waveforms are sent to a computer farm at the surface. In order to pass the trigger 
a minimum number of 8 DOMs in local coincidences within a time window of 5 //s is required. If 
this condition is fulfilled the waveforms are combined to an event and the number and arrival times 
of the Cherenkov photons are extracted. Here, the relative timing resolution of photons within an 
event is about 2 ns. The absolute time of an event is determined by a GPS clock to a precision of 
O(iis), which is more than sufficient for our analyses. 

Data are transferred from the South Pole to a computer center in the North via satellite. 
For the analyses described in this paper we consider only muons produced in charged current 
interactions 

v„ + N -> n + X (1) 

as only the track-like hit pattern of muons allows for a good angular resolution. In order to 
fit satellite bandwidth restrictions, a filter removes events which do not qualify as good upgoing 
neutrino candidates 4 . For this purpose, an initial track is reconstructed for each event using the 
line- fit algorithm (Ahrens et al. 2004). This is a simple but fast analytic track reconstruction based 
on the measured hit times in the DOMs. The transferred data forms the basis for our analyses. 

4. Data sets and reconstruction 

For our analyses we use data taken with the IceCube detector in its 22-string configuration 
from 2007 June 1 to 2008 April 4. To prevent bias in our analyses, the data within the — lh to 
+3h windows (on-time data) are initially only checked for detector stability until all parameters 
of the analysis have been fixed, i.e. they are not used in the optimization of the analyses. The 
remaining off-time data that pass basic quality criteria (95% of all the data collected in the 22-string 
configuration) amounts to 268.9 days of livetime. This long off-time window allows for a precise 



3 Six of these strings will make up a dense subarray in the clearest ice known as Deep Core, extending the sensitivity 
of IceCube at lower energies. 

4 A11 triggered data are stored on tape and shipped to the North during the Summer season. 
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experimental determination of the background rate in the on-time windows. 

In addition to the basic quality criteria, the on-time data is tested for stability by fitting a 
Gaussian to the distribution of low level event counts in each one second bin. In order to check for 
unexpected periods of dead time, an exponential is fit to the distribution of time delays between 
events. Data with deviations from the expected shapes are inspected more closely for possible 
causes which make the data unsuitable for our analysis. Here, most observed deviations can be 
attributed to transitions between runs and are therefore not critical. Out of 48 northern hemisphere 
bursts in the time period under investigation, the data for seven do not pass the stability/quality 
criteria or have gaps during the prompt /precursor emission windows. For all remaining 41 GRBs, 
both tests show excellent agreement with no indications of abnormal behavior of the detector during 
the on-time periods. The on-time data cover 100% of the prompt and precursor windows for the 41 
selected bursts. For the extended window seven out of the 41 bursts exhibit gaps at the beginning 
and/or end of the window. For these seven bursts the extended window is shortened accordingly 
(see Table 3). With this correction, the data cover 94% of the extended time windows of all GRBs. 
The missing 6% are due to larger gaps in data taking. 

After the data has been transferred to the North, a more precise determination of the direc- 
tion of an event is achieved by fitting a muon-track hypothesis to the hit pattern of the recorded 
Cherenkov light in the detector using a log- likelihood reconstruction method (Ahrens et al. 2004). 
A fit of a paraboloid to the region around the minimum in the log-likelihood function yields an 
estimate of the uncertainty on the reconstructed direction (Neunhoffer 2006). At this point, the 
data sample with an event rate of 3.3 Hz is still dominated by several orders of magnitude by misre- 
constructed downgoing atmospheric muons as demonstrated in Figure 2, which shows a comparison 
between data and Monte Carlo (see Section 5). The quantities shown are later used to reject these 
misreconstructed atmospheric muons and improve the sensitivity of the analyses. 

• 6> rec : reconstructed zenith angle 5 ; 

• Odir: the uncertainty on the reconstructed track direction (quadratic average of the minor and 
major axis of the la error ellipse); 

• L rc &: log of the likelihood value of the reconstructed track divided by the number of degrees 
of freedom (number of hit DOMs minus number of fit parameters). This has proven to be a 
powerful variable for separating signal and background as visible from Figure 2; 

• Lu/£,: difference in log-likelihood value between the reconstructed track and one containing a 
bias to be reconstructed as downgoing. The bias is zenith angle dependent and follows the rate 
of downgoing atmospheric muons. The rationale behind this is that a track is much more likely 
to originate from an atmospheric muon than from a muon generated in a neutrino interaction. 
Only high-quality upgoing tracks have high Ljj/ d values; 

• Adi r : the number of photons detected within a —15 to +75 ns time window with respect to the 



5 The zenith angle in detector coordinates is related to declination S by = 8 + 90°. 
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Fig. 2. — Comparison between data (black solid circles; 268.9 days of livetime) and simulations 
in the quality parameters used to reject misreconstructed atmospheric muons at filter level (see 
Section 3). Monte Carlo shown includes atmospheric muons (green solid lines), coincident muons 
(magenta dot-dashed lines), atmopheric neutrinos (red dashed lines), and prompt GRB neutrinos 
(blue dotted lines). The GRB signal is assumed to follow a standard Waxman-Bahcall spectrum 
and is normalized to the summed expectation from 41 bursts.. 
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expected arrival time for unscattered photons from the muon track hypothesis; 

• Omin. minimum zenith angle from a fit of a two-track hypothesis to the light pattern. For this, 
the light pattern is divided into two separate sets of hits based on the mean hit time of the 
event, and a track hypothesis is fitted to each hit set separately. A cut on the smaller zenith 
angle, 9 m m, of the two tracks is very effective against so-called coincident muons, where two 
muons from different atmospheric showers pass the detector in fast succession mimicking an 
upgoing track. In this case, often at least one of the two tracks is reconstructed as downgoing 
whereas for good-quality upgoing neutrinos both tracks appear most often as upgoing; and 

• E rec : reconstructed muon energy at point of closest approach to the center of gravity of hits in an 
event (Zornoza et al. 2008). This is the calibrated output of a likelihood reconstruction method 
evaluating the measured number of photons in each DOM with respect to the corresponding 
probability density function (PDF) of a given track-energy hypothesis. 

Muons carry a significant fraction of the original neutrino energy, and at the energies of highest 
acceptance in our analyses (~100TeV), have a range of about 10 km. The dominant energy 
loss mechanisms for these muons are Bremsstrahlung and pair production which grow with 
increasing energy, thereby increasing the amount of Cherenkov light emitted and allowing one 
to estimate the muon energy. 



5. Simulation 

Signal neutrinos are generated from the direction of the GRB with the corresponding spectrum 
using a port of the anis code (Gazizov & Kowalski 2005) called neutrino-generator. Here, the 
change in position of the source in the detector coordinate system during the respective time window 
is taken into account. The neutrino emission is assumed to be constant during this time. In addition, 
three types of background are simulated. At the beginning of the analysis chain, the data sample 
is dominated by downgoing atmospheric muons, produced by interactions of cosmic rays in the 
atmosphere, which are reconstructed as upgoing. As more cuts are applied the data sample starts 
to be dominated by coincident muon events. Both event classes are simulated with the CORSIKA 
air shower simulation package (Heck et al. 1998). Finally, we consider the irreducible background 
of atmospheric neutrinos generated in the same interactions as the atmospheric muons. These 
neutrinos are simulated as an all-sky flux with neutrino-generator and weighted according to 
the Bartol spectrum (Barr et al. 2004). 

Neutrinos are tracked from the surface through the Earth taking into account absorption, 
scattering and neutral current regeneration of neutrinos (Gazizov & Kowalski 2005). Information 
on the structure of the Earth is taken from the Preliminary Reference Earth Model (Dziewonski &: 
Anderson 1981). Muons originating from neutrino interactions near the detector and atmospheric 
muons are traced through rock and ice taking into account continuous and stochastic energy losses 
(Chirkin & Rhode 2004). The photon signal in the DOMs is determined from a detailed simulation 
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Table 3. Modified Extended Time Window for GRBs with Gaps in Data Coverage at Beginning 

or End of — 1 h to +3 h Window 



T[-T [s] n-T [ S ] 



GRB070610 


-1593 


+10800 


GRB070714B 


-3600 


+3930 


GRB071021 


-3600 


+5391 


GRB071109 


-3022 


+10800 


GRB080205 


-3600 


+6289 


GRB080211 


-3600 


+8399 


GRB080320 


-389 


+8440 



Table 4. Number of Signal and Background (Off-Time Data) Events for the Unbinned Method 

at Different Cut Levels 





Prompt Window 


Precursor Window 


Wide Window 


Cut Level 


No. Events 


e a (%) 


No. Events 


£ a (%) 


No. Events' 5 


e a (%) 




Signal 


Filter 


0.062 


100 


1.8 


100 


7.6 


100 


Final 


0.033 


53 


0.53 


29 


2.8 


37 


Background 


Filter c 


7.6 x 10 8 


100 


7.6 x 10 8 


100 


7.6 x 10 8 


100 


Final 


4846 


6.5 x 10~ 4 


4846 


6.5 x 10~ 4 


4846 


6.5 x 10~ 4 


Final (Window) d 


6.1 x 10~ 4 


8.2 x 10~ n 


6.1 x 10~ 4 


8.2 x 10~ n 


8.8 x 10~ 2 


1.2 x 10~ 8 



a EfBcicncy relative to filter level. 

b For a fluence equal to the computed upper limit in Section 7. 
C A11 events after cuts. 

d Expected events in cones with radii of 2.3° (contain 70% of signal events; see Figure 4) around GRBs 
within respective time window (prompt: T2 — Ti from Table 1; precursor: 100 s; wide: 4h). Note, that these 
cones are not used in the evaluation of the data with the unbinned likelihood method. 



-16- 



of the propagation of Cherenkov light from muons and showers through the ice (Lundberg et al. 
2007) which includes the modeling of the changes in absorption and scattering length with depth 
due to dust layers (Ackermann et al. 2006). This is followed by a simulation of the DOM electronics 
and the trigger. The simulated DOM signals are then processed in the same way as the data. 



6. Data analysis 

Muon neutrinos from GRBs show up as an excess of tracks above the background from the di- 
rection of the GRB within a certain time window. Background from misreconstructed atmospheric 
muons can be suppressed by applying quality cuts on reconstructed quantities. Background from 
atmospheric neutrinos on the other hand is indistinguishable from cosmic neutrinos. Hence, once 
a high-purity (atmospheric) neutrino sample has been selected, quality cuts which aim at reject- 
ing misreconstructed tracks cannot further improve the signal to background ratio. However, as 
neutrinos from GRBs are expected to exhibit a harder energy spectrum than that of atmospheric 
neutrinos, information on the muon energy allows to further increase the sensitivity of the analysis 
to neutrinos from GRBs. 

We have analyzed the data both with an unbinned likelihood and a binned method. In order 
to enhance the chances for a discovery, we do not analyze the 41 GRBs individually but as a 
population (stacked analysis) 6 . This allows us to set the most stringent limit on the tested models 
but might not be optimal in other cases, e.g. if one burst has a much higher neutrino fluence 
relative to the rest of the GRBs than expected. Also, in case of a discovery a stacked analysis only 
allows one to calculate the neutrino fluence from the whole burst population but not from individual 
GRBs. In the following we describe the unbinned likelihood and binned methods and compare their 
performances in the case of the prompt emission scenario. The more sensitive method is then used 
to obtain the results presented in Section 7. 



6.1. Unbinned log-likelihood method 

After filtering at the South Pole and transfer to the North the data sample is still dominated 
by downgoing muons which are reconstructed as upgoing. These muons are rejected by applying 
quality cuts on reconstructed quantities (for a description see Section 4). 



6 mc > 85° ; Odir < 3° ; 6> m in > 70° ; Lu/ D > 30 ; L rc d < < 



7.8 
8.5 
9.5 



for N dh < 7 
for iV dir = 7 
for N dir > 7 



(2) 



In case of the binned method, the numbers of expected and observed events of all GRBs are summed up. In case 
of the unbinned method, the stacking is performed in Equation (3). 
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After these cuts, a high-purity upgoing (atmospheric) neutrino sample remains with an event rate 
of 2.1 x 10~ 4 Hz. Table 4 lists the cut efficiency for data and signal. A comparison of data with 
simulations is displayed in Figure 3. In general, good agreement between the atmospheric neutrino 
Monte Carlo and data is observed. Small deviations for example visible at low N^h are most likely 
due to atmospheric muon background which is not accounted for by the Monte Carlo due to its 
limited statistics at this cut level. The cumulative point spread function is shown in Figure 4 (left). 
The median angular resolution is about 1.5°. Figure 4 (right) shows the muon neutrino effective 
area for different declination bands. 

The data sets after quality cuts are the starting point for the unbinned likelihood method. 
In contrast to binned methods where the event is rejected if it lies outside the cut region (binary 
selection), unbinned likelihood methods do not discard events but use PDFs to evaluate the prob- 
ability of an event belonging to signal or background population. The unbinned likelihood method 
used here is similar to that described in Braun et al. (2008). The signal, S(xi), and background, 
B(xi), PDFs are each the product of a time PDF, a directional PDF, and an energy PDF, where 
Xi denotes the directional, time, and energy variables. 

The directional signal PDF is a two-dimensional Gaussian distribution with the two widths 
being the major and minor axes of the la error ellipse of the paraboloid fit described in Section 4. 
The time PDF is flat over the respective time window and falls off on both sides with a Gaussian 
distribution. The width a of the Gaussian is determined by the length of the time window with a 
maximum of a = 25 s and a minimum of a = 2 s. The Gaussian accounts for possible small shifts 
in the neutrino emission time with respect to that of the 7-rays and prevents discontinuities in the 
likelihood function. The sensitivity of the method depends only weakly on the exact choice of a. 
The energy PDF is determined for each GRB individually. It is derived from the energy-estimator 
distribution of the tracks of the corresponding signal Monte Carlo data set (weighted to an E~ 2 
spectrum 7 ) after final cuts (see Equation (2)). The signal PDFs of the GRBs are combined using 
a weighted sum (Abbasi et al. 2006) 

Ef=i RBs Wj Sj(xi) 

Stot(Xi) = ^Ncm, ' ( 3 ) 

where Sj(xi) is the signal PDF of the jth GRB and Wj is a weight that in the case of the prompt 
and precursor window is proportional to the expected number of events in the detector according 
to the fluences described in Section 2. In the case of the extended window we use Wj = 1 for all 
GRBs in order to make the search as general as possible. 



7 Using an energy PDF different from that of the assumed signal spectrum cannot lead to an overestimation of the 
significance of a potential signal as the same PDF is also used to obtain the distribution of background-only samples 
(see Figure 5) from which the significance is calculated. Studies have shown that for our analyses the loss in discovery 
potential due to the different PDFs is very small. On the other hand, the softer energy PDF increases the sensitivity 
to scenarios where the true GRB spectra are softer than assumed. 
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Fig. 3. — Comparison between data (black solid circles; 15.99 days of livetime) and simulations 
in the quality parameters used to reject misreconstructed atmospheric muons at final cut level. 
The data includes all data-taking runs with an overlap with the (— lh to +3h) extended search 
windows (this may be considered the maximal on-time of the analysis). Monte Carlo shown includes 
coincident muons (magenta dot-dashed lines) , atmospheric neutrinos (red dashed lines) , and prompt 
GRB neutrinos (blue dotted lines). The simulated single atmospheric muons have been completely 
removed at this cut level and the statistics for the simulated coincident muons are very low. The 
GRB signal is assumed to follow a standard Waxman-Bahcall spectrum and is normalized to the 
summed contribution of 41 bursts. 
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Fig. 4. — Left: Cumulative point spread function for the unbinned method at final cut level. Right: 
Effective area for muon neutrinos in several declination bands as a function of energy after final 
event selection in the unbinned method. 



For the directional background PDF the detector asymmetries in zenith and azimuth must be 
taken into account. This is accomplished by evaluating the data in the detector coordinate system. 
The directional background PDF is hence derived from the distribution of all off-time events after 
final event selection in the zenith-azimuth plane of the detector. The time distribution of the 
background during a GRB can be assumed to be constant resulting in a flat time PDF. The energy 
PDF is determined in the same way as for the signal PDF with weights corresponding to the Bartol 
atmospheric neutrino flux. 

All PDFs are combined in an extended log-likelihood function (Barlow 1989) 

N 

\n(C((n s ))) = -(n s )-(n b ) + J2^(( n s)Stot(xi) + (n b )B(xi)) , (4) 

i=i 

where the sum runs over all reconstructed tracks in the final sample. The variable (725) is the 
expected mean number of background events, which is determined from the off-time data set. The 
mean number of signal events, (n s ), is a free parameter which is varied to maximize the expression 

in order to obtain the best estimate for the mean number of signal events, (n s ). 

To determine whether a given data set is compatible with the background-only hypothesis 10 8 
background data sets for the on-time windows are generated from off-time data by randomizing 
the track times while taking into account the downtime of the detector. For each of these data sets 
the ln(7£) value is calculated, yielding the distribution shown in Figure 5. The probability for a 
data set to be compatible with background is given by the fraction of background data sets with a 
larger \n(lZ) value. For comparison, the plot also displays the ln(lZ) distributions for background 
data sets with one and two injected Monte Carlo signal events, respectively. The signal events 
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Fig. 5. — Likelihood-ratio distribution of 10 8 randomized background-only data sets for the prompt- 
window analysis (solid line). Also shown are the corresponding distributions for background data 
sets with one (dotted) and two (fine dotted) signal events injected. The integrals of all distributions 
have been normalized to one. 

are randomly distributed among the GRBs, where the assignment probability to a specific GRB 
is proportional to the expected number of events from that burst. The energy of the neutrinos is 
generated according to the spectra calculated in Section 2.1. 



6.2. Binned method 

For the binned method, a machine learning algorithm was trained to separate signal and 
background. The algorithm used was a Support Vector Machine (SVM) (Cortes & Vapnik 1995) 
with a radial basis function kernel. It was provided with the best reconstructed track direction 
in detector coordinates as well as many quality parameters including o"di r , -^rcd, Lu/ D , N^ir, and 
m i n which are described in Section 4. The SVM was trained using the off-time filtered data as 
background and all-sky neutrino simulation weighted to the sum of the individual burst spectra 
as signal. The optimum SVM parameters (kernel parameter, cost factor, margin) were determined 
using a coarse, and then fine, grid search with a 5-fold cross validation technique at each node, as 
described in Hsu et al. (2003). 

The resulting SVM classification of events is shown in Figure 6. The final cut on this classifier 
is optimized to detect a signal fluence with at least 5a (significance) in 50% of cases (power) by 
minimizing the Model Discovery Factor (MDF) according to Hill et al. (2006). The MDF is the 
ratio between the signal fluence required for a detection with the specified significance and power 
and the predicted fluence. The angular cut around each GRB is then calculated to keep 3/4 of 
the remaining signal after the cut on the SVM classifier. In this way, there is one cut on the 
SVM classifier for all GRBs, but different angular cuts around each GRB according to the angular 
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Fig. 6. — The SVM classifier distribution of data (black solid circles) and simulations. Monte Carlo 
shown includes atmospheric muons (green solid lines), coincident muons (magenta dot-dashed lines), 
atmospheric neutrinos (red dashed lines), and prompt GRB neutrinos (blue dotted lines). The GRB 
signal is assumed to follow the summed calculated individual neutrino spectra and is normalized 
to the rate of atmospheric neutrinos. The vertical dashed line indicates the final optimum cut at 
0.25. 

resolution of the detector in that direction. The optimum SVM cut is determined to be at a value 
of 0.25. Table 5 displays the signal and background event rates at different cut levels. 




6.3. Comparison of the two methods 

We compare the performance of the unbinned likelihood and binned method by means of their 
discovery potential for the prompt neutrino emission scenario. Figure 7 displays for a significance 
of 5(7 the power as a function of the MDF for the two methods. For a power of 50% the unbinned 
method shows an improvement in the MDF of about a factor 1.8 compared to the binned method. It 
is therefore used to derive the results presented in the following section. The large gain in sensitivity 

Table 5. Number of Signal and Background (Off-Time Data) Events for the Binned Method in 

the Prompt Window at Different Cut Levels 







Signal 


Back 


ground 


Cut Level 


No. Events 


Efficiency 3 (%) 


No. Events 


Efficiency 3 (%) 


Filter 


0.062 


100 


77 x 10 7 


100 


Final 


0.023 


37 


4.7 


6.1 X 10~ 9 



a Rclativc to filter level. 
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Fig. 7. — Comparison of the discovery potentials of the unbinned likelihood (solid) and binned 
(dashed) method for significances of 3.3a (light) and 5<r (dark). Shown is the fraction of data sets 
yielding at least the stated significance (power) as a function of the ratio between the mean number 
of injected signal events and the expected number of signal events from the model (MDF). 

is partly due to the explicit use of energy information in the likelihood. The discovery potential for 
the expected fluxes is further improved by weighting the bursts according to the expected number 
of signal events in the detector in the unbinned method (Equation (3)). In the binned method all 
bursts are treated equally. 
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7. Results and systematic uncertainties 

We apply the unbinned likelihood method to the on-time data sets after neutrino candidate 
event selection with the final cuts (Equation (2)). For all three emission scenarios the values of 
\n(lZ) and (n s ) are zero and hence consistent with the null hypothesis. Therefore, we derive 90% CL 
upper limits 8 on the fluence from the 41 GRBs in the prompt phase of 3.7 x 10 -3 erg cm -2 (72 TeV 
- 6.5 PeV) and on the fluence from the precursor phase of 2.3 x 10 -3 erg cm -2 (2.2 TeV - 55 TeV), 
where the quoted energy ranges contain 90% of the expected signal events in the detector. Further 
information is listed in Table 6. The limits, which are displayed in Figure 8, are not strong enough 
to constrain the models. The 90% CL upper limit for the wide time window is 2.7 x 10 -3 erg cm -2 
(3 TeV - 2.8 PeV) assuming an E~ 2 flux. 

To illustrate, we counted the number of events after final cuts in cones with radii 2.3° around 
the GRB positions (contain 70% of signal events; see Figure 4) within the corresponding time 



8 The limits are calculated following a procedure due to Neyman (Neyman 1937; Amsler et al. 2008). The 90% CL 
upper limit corresponds to the signal flux for which 90% of background data sets with signal events injected according 
to the calculated GRB spectra yield ln(72.) values greater than the observed one. 
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Fig. 8. — 90% CL upper limits (dark thick lines) on the neutrino fluence from the 41 northern 
hemisphere GRBs for different emission models (light thick lines): precursor (dashed, Razzaque 
et al. (2003)) and prompt (solid, see Section 2). The dark dotted and dash-dotted thin lines mark 
the scaled AMANDA 90% CL upper limits (Achterberg et al. 2008; Kuehn 2007) on the prompt 
and precursor nuences, respectively. 

windows (prompt: T<i — T\ from Table 1 ; precursor: 100s; wide: 4h). In addition, we analyzed 
the data in the prompt window with the binned method after final cuts. In all cases, zero events 
remain which is consistent with the results of the unbinned likelihood method. 

As described previously, we use the off-time data to determine the background rate in the on- 
time windows. This technique removes many potential sources of uncertainties in the calculation of 
the significance of a possible signal that are introduced when using a simulation of the background. 
However, this method makes the assumption that the rate of data during the off-time and on- 
time windows are the same. Furthermore, we use Monte Carlo for the signal simulation and the 
derivation of upper limits, which involves the propagation of particles through the Earth and ice, and 
the simulation of the detector response. The most important sources for systematic uncertainties 
are discussed below in detail. Their effects on the upper limits are summarized in Table 7. 

Ice simulation: Inaccuracies in the ice simulation can lead to a wrong estimate of the efficiency 
of the detector to neutrinos from GRBs. Data-Monte Carlo comparisons and variation of simulation 
parameters indicate that the systematic uncertainty from this aspect of the simulation is about 



±15%. 
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Table 6. Summary of Search Results for Prompt and Precursor Window 



Window 


^-exp 


"limit 


Factor 


Prompt Window 


0.033 


2.4 


72 


Precursor Window 


0.26 


2.5 


9.7 


Wide Time Window 




2.7 





Note. ~~ i-cxp : number of expected events in 
the detector after final cuts from all 41 GRBs in 
the unbinned search; nii m i t : 90% CL upper limit 
on the event number from all 41 GRBs; factor: 
factor by which the limit exceeds the predicted 
event number n exD . 



Table 7. Summary of Effects of Systematic Uncertainties on the Upper Limits 



Type of Uncertainty 


Prompt Window 


Precursor Window 


Extended Window 


Ice Simulation 


±15% 


±15% 


±15% 


DOM Efficiency 


±5% 


±10% 


±7% 


Lepton Propagation 


±5% 


±5% 


±5% 


Background Rate 


< 1% 


< 1% 


< 1% 


Sum 


±17% 


±19% 


±17% 
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DOM efficiency: A ±10% uncertainty in the efficiency of the optical modules in the detection 
of photons leads to a corresponding uncertainty in the number of expected events from a GRB. 
This effect is nonlinear and spectrally dependent, so simulation was generated spanning the range 
of uncertainties to determine the resulting change in signal event rates; 

Neutrino and muon propagation: Theoretical uncertainties on muon energy losses and the 
neutrino-nucleon cross-section, determined from the uncertainty on the CTEQ6 PDFs (Pumplin 
et al. 2002), contribute a 5% uncertainty on the neutrino event rate in the detector; 

Background rate: After final cuts the variation of the event rate over the data taking period is 
about ±5%. In order to account for potential differences at the time of the bursts the background 
data rate is varied by this amount. This results in a shift of the upper limits of less than ±1% and 
is therefore negligible. 

8. Comparison to other results 

A search with the AMANDA detector for muon neutrinos in the prompt and precursor phase 
was conducted for GRBs detected between 1997 and 2003 (Achterberg et al. 2008) with null re- 
sult. The analysis of prompt neutrinos contained 419 bursts observed by the BATSE experiment 
(Paciesas et al. 1999) as well as by others with similar characteristics. The large number of bursts 
allowed that analysis to set an upper limit only a factor 1.4 above the prediction of the Waxman- 
Bahcall prompt emission model. A further search with the AMANDA detector for muon neutrinos 
was conducted for 85 GRBs detected between 2005 and 2006 (Strahler 2009), primarily by the 
Swift satellite, also with null result. Due to the smaller number of bursts, the upper limit from 
this analysis is much less restrictive. Converting the upper limit obtained from the 419 bursts to 
a fluence limit from 41 standard Waxman-Bahcall bursts (for definition see footnote 2) yields the 
dotted line in Figure 8. Due to the 10 times smaller number of bursts available, the limit presented 
in this paper is about a factor three worse. For the precursor emission model, the AMANDA limit 
(Achterberg et al. 2008; Kuehn 2007) is much less restrictive as only 60 bursts, detected between 
2001 and 2003, were used. It is shown in Figure 8 as a dash-dotted line. Here, our analysis improves 
on the AMANDA upper limit despite the fact that 30% fewer bursts were investigated. 

The AMANDA data were also analyzed for neutrinos of all flavors from GRBs (Achterberg 
et al. 2007). Apart from a search for neutrinos from 73 bursts detected by BATSE in 2000, another 
search did not rely on information from satellites but looked for a clustering of events within sliding 
time windows of 1 s and 100 s. Due to this more generic approach and the low number of bursts, 
respectively, the limits from these analyses are much less restrictive than those from Achterberg 
et al. (2008). However, the sliding-window or similar searches are the only way to detect GRBs 
where either the jet does not emerge from the progenitor star (choked bursts, Meszaros & Waxman 
(2001)) or the 7-ray signal is not observed by satellites for other reasons. 

We continue with a discussion of the impact of results from high-energy 7-ray observations on 
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expected neutrino fluences. Along with high-energy neutrinos which originate from the decay of 
charged pions, high-energy 7-ray photons are produced in the decay of simultaneously generated 
neutral pions. In addition, high-energy photons are produced in inverse-Compton scattering of 
synchrotron photons by accelerated electrons (Falcone et al. 2008). In contrast to neutrinos, the 
flux of high-energy 7-ray photons at the Earth is significantly reduced due to the large optical 
depths for photon-photon pair production inside the source for not too large jet Lorentz factors 
Tjct < 800 (Falcone et al. 2008). In addition, high-energy photons above 100 GeV are absorbed on 
the extragalactic background light (EBL) if they travel distances with z > 0.5. Observations with 
air-Cherenkov telescopes like H.E.S.S. (Aharonian et al. 2009; Tarn et al. 2009) or MAGIC (Albert 
et al. 2007; Galante et al. 2009) are also hampered by the fact that usually it takes more than 
50 s (MAGIC) or 100 s (H.E.S.S.) from the observation of a GRB by a satellite to the start of data 
taking with these telescopes. Therefore, the prompt emission window is only partially covered or 
not at all. MILAGRO as an air shower array observed large parts of the sky continuously (Atkins 
et al. 2004, 2005; Abdo et al. 2007). However, it was mostly sensitive to energies above 100 GeV 
and therefore suffered significantly from 7-ray absorption on the EBL. HAWC, the successor of 
Milagro currently in the planning phase, will be able to detect 7 rays from GRBs down to 100 GeV 
where about 50 (Fermi) GRBs per year will fall in its field of view (Goodman et al. 2009). 

At energies above 100 GeV, there has been no definitive detection of 7-ray emission from GRBs. 
Milagrito (Atkins et al. 2000, 2003) and the HEGRA AIROBICC array (Padilla et al. 1998) reported 
evidence at the 3a level for high-energy 7-ray emission from GRB 970417A (E 1 > 650 GeV) and 
GRB920925C (E^ > 20TeV), respectively. However, subsequent searches for high-energy 7-ray 
emission from GRBs did not find similar signals. The limits obtained from MAGIC and H.E.S.S. 
are not directly comparable to our results due to their incomplete burst coverage. In Atkins et al. 
(2004), Atkins et al. (2005) and Abdo et al. (2007) the MILAGRO collaboration reports 99% 
CL upper limits for a large number of individual bursts (both long and short) between 2000 and 
2006 down to 10~ 7 erg cm -2 (energy range ~100 GeV - 10 TeV; the exact energy range differs from 
publication to publication). Extrapolating the average per burst 99% CL upper limit for the prompt 
window in our analysis to the energy range from 100 GeV to 10 TeV yields 3.2 x 10 -7 ergcm~ 2 . 
However, the photon limits do not constrain our results as they do not account for absorption in 
the EBL (this would significantly worsen the limits) and include an unknown, probably dominant, 
contribution from inverse-Compton scattering. In general, current flux predictions for high-energy 
gamma rays from GRBs are near or below the sensitivity of current instruments (Falcone et al. 
2008), where the predicted fluxes in the energy range below ~ 100 TeV are dominated by the leptonic 
emission component in most scenarios. 

Within the internal shock (fireball) model, synchrotron self-Compton (SSC) processes between 
the accelerated electrons and the 7-ray photons could lead to a high-energy 7-ray peak in the 
GeV range detectable by Fermi in case of bright GRBs. Up to now, this has happend only for 
a handful of bursts (one of the photons with the highest energy from a GRB detected by Fermi 
came from GRB080916C and had an energy of ~ 13 GeV (Abdo et al. 2009)). However, this does 
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not necessarily disfavor the internal shock scenario. As discussed in Fan (2009), the amount of 
energy in high-energy photons could be suppressed by an inefficient SSC process in the extreme 
Klein-Nishina regime or a combination of a SSC peak at high energies and a low photon cut-off 
energy above which the fireball becomes optical thick. 

9. Conclusions and outlook 

We have performed a set of complementary searches for muon neutrinos associated in space 
and time with 41 gamma-ray bursts that were observed in the northern sky between 2007 June 
and 2008 April. For the first time in searches with large GRB populations, we have calculated 
individual prompt neutrino spectra for all 41 GRBs using measured GRB parameters. The search 
results are consistent with the case of a background-only hypothesis. Therefore, we place 90% CL 
upper limits on the fluence from the prompt phase of 3.7 x 10~ 3 erg cm -2 (72 TeV - 6.5 PeV) and on 
the fluence from the precursor phase of 2.3 x 10~ 3 ergcm~ 2 (2.2 TeV - 55 TeV), where the quoted 
energy ranges contain 90% of the expected signal events in the detector. Though the number of 
bursts is smaller than in previous searches the larger detector allows us to improve on the limits for 
the precursor phase by a factor 1.4. Compared to the predictions, the limits lie a factor 72 (prompt 
phase) and 9.7 (precursor phase) higher. Hence, they do not allow us to constrain the models. 
Apart from these model-driven searches, we have also conducted for the first time a generic search 
for neutrino emission from GRBs in a wide window of (— lh to +3h) around each burst. Finding 
no evidence for a signal, we place a 90% CL upper limit on the fluence of 2.7 x 10~ 3 erg cm" 2 (3 TeV 
- 2.8 PeV) assuming an E~ 2 flux. 

Launched in 2008 June, the Fermi Gamma-ray Space Telescope (NASA 2008) has begun to 
provide an expanded catalog of sources for future neutrino searches. With a much larger field of 
view than other satellites, Fermi has increased the GRB detection rate by more than a factor two. 
At the same time, the detected bursts have on average a higher luminosity than those detected 
by Swift due to the lower sensitivity of the Fermi-GBM (Gamma-ray Burst Monitor) instrument. 
The first IceCube analysis to take advantage of the increased detection opportunities will utilize 
the 40-string configuration of the detector, already as large as the full IceCube along one axis. 
This gives it the full angular resolution power along that direction and thus provides powerful 
background rejection. However, this is mitigated by the comparatively poor angular resolution 
of the Fermi-GBM (~ 3°), which is worse than the IceCube resolution. With the full 80-string 
detector scheduled to be completed in 2011 and an expected 100 to 150 detected bursts per year in 
the northern hemisphere by the Fermi and Swift satellites the sensitivity of IceCube to neutrinos 
from GRBs will soon exceed that of AMANDA. This will allow IceCube to either confirm the 
predicted fluxes within the next years or set stringent limits thereby disfavoring GRBs as the 
major sources of ultra-high energy cosmic rays. 
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• Parameters of the 7-ray spectrum F 1 {E 1 ): 

- e 7 : break energy; 

- q 7 : spectrum index before break energy; 

- /3 7 : spectrum index after break energy; 

- Ty. measured fiuence in 7-rays integrated over the energy range given in the GCN circulars 
and reports (NASA 2009); 

- / 7 : normalization; obtained from integral of Equation (A2). 

• Parameters of the neutrino spectrum F U (E U ): 

- €\\ first break energy; 

- 62- second break energy; 

- a u : spectrum index before first break energy; 

- (3 V : spectrum index between frist and second break energy; 

- 7j,: spectrum index after second break energy; 

- f v : normalization; obtained from integral of Equation (A8). 

• z: redshift of GRB; 

• e e : fraction of jet energy in electrons; 

• €b- fraction of jet energy in magnetic field; 

• f e : ratio between energy in electrons and protons; 

• L™°: isotropic luminosity of the GRB; 

• t var : variability of the 7-ray light curve of the GRB; 

• Tjct: Lorentz boost factor of the jet. 

The expression 1 — (1 — (x p _ +7r )) A ^ Ap7 in Equation (A8) estimates the overall fraction of the 
proton energy going into pions from the size of the shock, AR, and the mean free path of a proton 
for photomeson interactions, A P7 . Here, (x p -+ n ) = 0.2 is the average fraction of proton energy 
transferred to a pion in a single interaction. The expression ensures that the transferred energy 
fraction is < 1. The calculations are insensitive to the beaming effect caused by a narrow opening 
angle of the jet as all formulae contain the isotropic luminosity in conjunction with a Air shell 
geometry, i.e. effectively use luminosity per steradian. For example, the target photon density used 
to calculate N- Uit is given by n 7 oc L 7 ° /4irR 2 , where R is the distance of the shock region from the 
central black hole. 
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